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ABSTRACT 



The use of Laplace Transforms has long been one of the primary 
methods utilized in solving linear differential equations but recently the 
extension and application of this method to the solution of basic non- 
linear systems has been proposed. This latter technique is presented 
and expanded to include various non-linear functions . The applicability 
of this method to systems of varied form and complexity is then explored 
and the validity of the derived solution investigated. 
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I. INTRODUCTION 



Laplace Transforms have been typically applicable and limited to 
linear systems. However, Baycura [Ref. 1 and 2 ] in two recent papers 
has proposed an extension of the Laplace Transform method to non-linear 
systems. By expanding the Laplace integral into an infinite series and 
integrating by parts n - times an infinite series was developed for non- 
linear terms. Through a further correlation of these derived series with 
series of linear transforms the transform for a non-linear term was obtain- 
ed. These papers develop this technique for higher-order derivatives and 
powers of the dependent variables illustrated with several applicable 
examples . 

The purpose of this paper will be to present the general method as 
developed by Baycura and to extend this method to varied and more com- 
plicated equations of non-linear systems . The latter will be accomplished 

by developing the transforms for the product of derivatives of the depen- 

th 

dent variable and the dependent variable raised to the n power, and 
derivatives raised to the n^ 1 power. The validity of the transforms will 
then be shown in several applicable examples. In the examples a 
computer method designed to find the inverse Laplace Transform will be 



introduced and explained. 



II. LAPLACE TRANSFORM CONCEPT 



A. THE METHOD OF INTEGRATING BY PARTS 



The following generally summarizes Baycura's method of utilizing 
Laplace Transforms on non-linear systems using the method of integration 
by parts. Considering the product of two time-dependent functions, x(t) 
and y(t) , which are to be integrated using integration by parts n - times, 
obtain the general formula: 



jo uu go uo 

xydt = x J ydt - & J ydtdt + X J'J' ydtdtdt - . 

) 0 o o 

00 

(-l)n - 1 x (n-i) J ( J ydt n-i^ dt + . . . . 



. • + 



( 2 - 1 ) 



n = 1 , 2 , 3 , 



where x^ n_1 ^ indicates the (n - l)^ 1 derivative of x. 

Applying this general formula to the Laplace Integral, equation (2-1) 

can be used to find the series form of the Laplace Transform of a function 

x(t) allowing y(t) = e -st . This integration gives the series form of: 

00 00 00 00 



xe 



-St 






>-st 



dt - x 



J'e st dtdt+M Jj J e st dtdtdt-... (2-2) 



© 



0 



Evaluating at the upper limit it is found that the value of the integral 

approaches zero as t approaches infinity. Then evaluating at the lower 

limit the integral becomes: 

00 



J xe 



-St 



m 

s 



+ 



m 



+ 



m 

„3 



(2-3) 



0 



or 
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£[x(t)] = x(s) = + . 



w • • 



(2-4) 



B. APPLICATION TO DERIVATIVES 

Equation (2-4) depicts the Laplace Transform of x(t) in an infinite 
series expansion involving the initial conditions of the function. Through 
similiar manipulations the transforms of the derivatives may be obtained: 




= sX(s) - x(0) (2-5) 



00 

£ [5i(t)3 = jx e" st dt 



0 



m ) 

s 



+ 




+ 



= s 2 X(s) - sx(0) - x(0) 



(2-6) 



C. APPLICATION TO NON-LINEAR FUNCTIONS 

The simplest non-linear function to investigate is x 3 (t) . Using the 
same technique, expansion by parts yields: 



£ [x 2 (t) ] = x x e s " dt - * J x e st ' dt dt + X 



J'J 



xe st dt dt dt 



(2-7) 



Taking each of the above terms individually and again expanding by parts 
obtain: 



xjxe st dt = x(0) X(s) 



- x 



J xe- st dtdt = ^LXfet 
J s 



( 2 - 8 ) 



(2-9) 
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X 



J J x e -st dt dt dt = 



_ X( 0) X(s) 



( 2 - 10 ) 



The summation of these terms then yields: 



£[x 3 (t )] = [x(0) + 



mi ^ mi 



+ ] X(s) = sX 3 (s) (2-11) 



Through similiar computations , transforms of higher-order functions 
may be obtained. The transform of the function cubed would be: 

X[x 3 (t)] = [x(0) + ^ ^ + ... ] s X(s) = s 3 X 3 (s) (2-12) 

Generalizing it can be stated that: 



X [x n (t)] = s n_1 X n (s) 



(2-13) 
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III. LAPLACE TRANSFORMS FOR NON-LINEAR FUNCTIONS 



A. TRANSFORMS FOR THE PRODUCT OF A FUNCTION AND ITS FIRST 
DERIVATIVE 

Transforming x(t) x(t) by applying integration by parts the following 
derivation is obtained: 



X [*(t)x(t)] = 



xxe st dt = x 



x e st dt - 3c 



J x e -st dt dt 



(3-1) 



+ x 



JJ’xe st dtdtdt - ... 



Let: I x = x j x e st dt 



Knowing the transform of the term under the integral, this term then yields 
I x = x(0) X(s) (3-2) 

Similiarly allow: 



Ip r -xJjxe st dtdt = x(0) 



X(s) 



(3-3) 



and: 



In = x 



JJxe st dtdtdt= x(0) 



Xls) 



(3-4) 



The summation of these terms gives: 



X [xx] = E I n = [Ii + I 8 + I 3 + . .. ] 
n=I 



= x(s) [=5(0) + + 



+ . 



(3-5) 
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Simplifying: 



io«Dx(t)]-s'xw[-sjia + + ^ + ^ + ...] 



= s 2 



X(s) jx(s) - f* ] 



= sX(s) [s X(s) - x(0)] (3-6) 

Continuing in this fashion the Laplace Transform for x(t)x a (t) is 



derived. 



£[x(t) x 2 (t)] = 



^ xx 2 e st dt = x 



00 co 

x s e -st dt - k Jj x a e _st dtdt 



+ x 



l ’x 3 e st dtdtdt - ... 



(3-7) 



Let: I x = kjx 2 e st dt = xjxjxe -st dt - xjj'xe st dtdt + .... 



Further let: J x = x x e -st dt = x(0) X(s) 



I, - 



Ja * x jj'j' xe" st dtdtdt = 



Summing: 



Ii = k I =1 J n = *(0) X(s) [x(0) + ~“ + 
and after simplifying: 



MO) j. MQ] 

s 2 
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I x = s x(0) X s (s) 



(3-8) 



Similiarly let: 



- _ V PP ^ 2 ^-St 



- - 



x e -bl dtdt = 






xe st dtdt - x 



JJ xe~ st dtdtdt + . . .J 



L-, = x 



xe” st dtdt = - 
J s 



L 3 =-x j|rxe- st dtdtdt»-^5^ 



Summing again: 

= x(0) X 2 (s) 



(3-9) 



By similiar arguments obtain: 



1.9 ~ 



- 2 j If 



x 3 e st dtdtdt = 



_ x(0) X 3 (s) 



(3-10) 



Combining the above integral terms the transform is then: 

x tx(t) x 3 (t)] - £ . I n = s 3 X s (s) r. 210 J + sM + + . 

n - 1 l s s s 



•] 



= s 2 X 2 (s) [s X(s) - x(0) ] (3-11) 

In similiar fashion the transform for x(t) x 3 (t) may be obtained and 
expressed as : 



00 

£ [x(t) x 3 (t)] = Jxx 3 e St dt = s 3 X 3 (s) [s X(s) -x(0)] 



(3-12) 



Equations (3-6), (3-11) and (3-12) provide the Laplace Transform for 
the general form of the function, x(t) x n (t). 
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£[*(t) x n (t) J = s n X n (s) U[*(t)] } 



(3-13) 



B. TRANSFORMS FOR THE PRODUCT OF A FUNCTION AND ITS SECOND 
DERIVATIVE 

Employing again the method of integration by parts to the non-linear 
function, x(t)x(t), the transform may be expressed as: 



£[x(t)x(t)] = Jxxe -St dt = xjxe' 



st 



dt - x 



J xe -st dtdt + ...(3-14) 



Attacking each of the integrals separately as in Section III. A we let: 

00 

I x = xjxe -st dt = X(0) X(s) (3-15) 



Ip. = - x 



f Jxe -st dtdt = aoms) 

0 c 



la = X v 



Therefore: 



= vW ijjxe _st dtdtdt = X - M * is) 



( 4 ), 



(3-16) 



(3-17) 



X [x(t)x(t)0 = 2 =1 I n = x(0) X(s) + + s (4) (0^ X(s) + ... 



= s 



X(s) [- 



_ x(0) _ x(0] x(0) x(0) X(0} 



. . .] (3-18) 



5 S S S S 

Collecting terms and simplifying: 

£[X(t) x(t)] = s X(s) [s a X(s) - sx(0) - £(0)] (3-19) 

The derivation of the Laplace Transform of X(t)x 2 (t) may be similiarly 
obtained . 



fx x 3 e dt = X 


>e" st dt - x f 


J J 


J' 



Jx 2 e st dtdt + ... (3-20) 
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Knowing the transform of x 3 (t) again write separately for each of the 
integrals : 



Ii = X 



x 3 e -st dt = k’(0) sX 2 (s) 



(3 



la = " X 



J'x 2 e _St dtdt = x(0) X 2 (s) 



(3 



I 3 = x( 4 ) pJx B e st dt dt dt = 



(3 



Summing the above terms the transform becomes: 

i[X(t)x 2 (t)] = | =1 I n = X(0)sX 3 (s) + x(0) X 3 (s) + sbMj£kl + .. 

= s‘X=(s)r-^-^ + ^ + ^f + ^ + ...l (3 

S S S S S _) 



£[k(t) x 2 (t)]= s 2 X 2 (s) [s 2 X(s) - sx(0) - x(0)] 

The transform for k(t)x 3 (t) may be identically expressed as: 



(3 



00 

r» 



£[k*(t) x 3 (t)] = j xx 3 e -st dt =k 



r x 3 e st dt - x 



Jx 3 e -St dtdt + . . . (3 



£[k(t) x 3 (t)] = s 3 X 3 (s) [s 2 X(s) -sx(0) - *(0)] (3 

Comparison of equations (3-19), (3-25) and (3-27) provides the 
general form of the transform for these functions. 



£[k(t) x n (t)] = s n X n (s) [s 2 X(s) - s x(0) - *(0)] 
= s n X n (s) {£ [x(t)]} 



(3 



- 21 ) 

- 22 ) 

-23) 

-24) 

-25) 

-26) 

-27) 

-28) 
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c. TRANSFORM OF THE PRODUCT OF POWERS OF A FUNCTION AND ITS 
DERIVATIVES 

Combining the results of the two foregoing sections the general form 
of the transform of the product of a function raised to the n^ 1 power and 
its various derivatives can be derived. From the derived transforms as 
expressed in equations (3-13) and (3-28): 

I[*(t) x n (t) 3 = s n X n (s) {£ [X(t)]} 

X[x(t)x n (t)] = s n X n (s) UDc(t)]} 

The Laplace Transform of the product of a function raised to the n^ power 
and its m 1 -* 1 derivative is therefore: 

I[x( m )(t)x n (t)] = s n X n (s) [£[x (m) (t)]} (3-29) 

D. THE TRANSFORM OF THE FIRST DERIVATIVE RAISED TO THE n th POWER 
From Section II-B, equation (2-5) gives the transform of x(t). 

£ [&(t) ] = s X(s) - x(0) 

Proceeding with the first derivative squared, [x(t)] 2 , and expressing 
the transform using integration by parts, the Laplace Transform becomes: 



£[(x) 2 ] = 



x e 



~ st = X f 






x e -st dt 



J - Xe -St dtdt + 



(3-30) 



0 



0 



Knowing the transform of X(t) the terms of the series may be expressed 
separately. 

00 



h = X j"x e -st dt = x(0) [s X(s) - x(0)] 



(3-31) 



0 



Is = -X 



00 

J X e -st dt dt = [sX(s)-x(0)] 



(3-32) 
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I a = X 



"JJ xe _st dtdtdt = [sX(s) -x(0)] 



_ x(0) 



(3-33) 



£ { [x(t)] 2 } = | =1 I n = I x + I a + I a + . . . 



= [s X(s) - x(O)] jx(0) 



+ m + m + 

s s 2 






= [sX(s)-x»=[-f 
- s s CsX(s) - x(0)] [xls)- 2 ^; 

= s[sX(s) - x(0)] 2 (3-34) 

Using similiar arguments the transform of [£(t)] 3 may be derived. 
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£ {x(t)] 3 } = J& 3 e -St dt = xjx 3 e -st dt -X Jx 2 e st dtdt+ ... (3-35) 

0 0 9 

Knowing the transform of [x(t)] 2 from equation (3-34) let: 



I x = x J* 2 e st dt = x(0) s [s X(s) - x(0) ] s 



(3-36) 



1 2 = - 5cjj:& 2 e -st dtdt = s [s X(s) - x(0)] 2 

0 

06 #t , 

1 3 = x J]‘x 2 e -St dtdt dt = s [sX(s) - x(0)] : 

J s 

0 

Summing the above terms: 

£ {[x(t)] 3 }= £ =1 I x + I s + I 3 + ... 



(3-37) 



(3-38) 



= s [s X(s) - x(0)] 2 j*(0) + + . . . I 

L s s J 
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= s 3 [s X(s) - x(0)] 2 [- 



x(0] + x(0) + MO] + x^O] + 
s s s 2 s' 3 



= s 2 [sX(s) - x(O)] 3 (3-39) 

From equations (3-34) and (3-39) the transform of [x(t)] n may be obtained. 

£ {[x(t)] n } = s 11-1 [s X(s) - x(0)] n 

= s"- 1 {X[x(t)]} n (3-40) 
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IV. APPLICATION OF LAPLACE TRANSFORMS TO NON-LINEAR SYSTEMS 



A, EXAMPLE OF A FIRST-ORDER NON-LINEAR SYSTEM, EQUATION 1 

The first non-linear system to which the derived Laplace Transforms 
are to be applied is one which has previously been considered by Baycura 
[Ref. 3] and Brady [Ref. 4] in the initial work that has been accomplished 
in this area. The differential equation to be considered is of the form: 

x + ax + bx 2 = 0 

The reasons for using this particular equation, aside from providing 
some measure of continuity in previous work, are that it represents a 
stable system and is by far one of the simplest non-linear equations 
with which to work. 

In this initial analysis the author will bring to light several concepts 
and techniques that utilize Laplace Transforms to derive solutions to 
various equations and then compare these to solutions obtained from more 
"conventional" techniques. The general concept will be to express the 
transform of the unknown variable in a series of terms in s , the solution 
to which may easily be obtained on the digital computer. The two solu- 
tions used for comparison will have been obtained from an analytical 
method and by computer using Runge-Kutta fourth-order integration methods 
using Adams -Moulton Predictor Corrector with error check. 

Utilizing, therefore, the equation: 

x + ax + bx 2 = 0 (4-1) 

and applying the derived Laplace Transforms we obtain: 
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bsX 2 (s) + (s + a)X(s) -x(0) = 0 



(4-2) 



Letting x(0) = c then: 

bsX 2 (s) + (s + a)X(s) - c = 0 
Using the quadratic formula, solve for X(s). 



(4-3) 



X(s) = 



s + a 



{-‘±[ 



1 + 



4bcs 



l 



] } 



2bs l — L/ ' (s+a) 3 

The radical may be expanded in a binomial series of the form of: 



(1+M) n = lt„M + ... 



(4-4) 



(4-5) 



Letting a = 2, b = 1, and x(0) = c = 1 the series expression for X(s) is: 



X(s) = 



1 



s+2 (s+2) 3 (s+2) 5 (s+2) 7 (s+2) e (s+2) 



2s‘ 



5s‘ 



14s 4 



42s & 132s 

u + 



(s+2) 13 •' 

(4-6) 



To obtain x(t) the inverse transform of the individual terms may be obtain 
ed and summed. An alternate and by far a more convenient and simpler 
method was, however, used in which X(s) was expressed as a ratio of 
polynomials in s. X(s) then took the form: 



X( s ) = bj^ 11 " 1 +b 2 s n ~ 2 + .... +b n 

s n + aiS 11 ”! + a 2 s n “ z +. . . . + a n 

A devised computer technique was then used to obtain x(t). This method 
to find the inverse Laplace Transform by digital computer is thoroughly 
described in Appendix A. 

Using these computer methods a solution to the equation was easily 
obtained using any desired number of terms of the series. As a result, 
solutions to the equation were examined terminating with the second to 
the eighth term in the series . 

18 



Several sets of solutions have therefore been derived and are tabulat- 
ed in Table IV- 1. These are an exact solution which is known and obtain- 
ed through analytical methods [Ref. 5], a computer solution using Runga- 
Kutta integration methods , and several computer solutions using Laplace 
Transforms but using varied numbers of terms in the truncated series. The 
Runga-Kutta solution is presented in order to show its accuracy since this 
solution will be the basis for comparison in subsequent examples. 



TABLE IV -1 

SOLUTION TO EQUATION 1 

EXACT RUNGE-KUTTA APPROX. APPROX. APPROX. 



TIME(t) 


SOLUTION 


SOLUTION 


4 TERMS 


6 TERMS 


8 TERMS 


0.0 


1.000 


1.000 


1.000 


1.000 


1.000 


0.2 


0.576 


0.575 


0.579 


0.579 


0.582 


0.4 


0.352 


0.352 


0.367 


0.369 


0.370 


O 

CD 


0.223 


0.223 


0.247 


0.249 


0.252 


1.0 


0.0950 


0.0944 


0.131 


0.101 


0.140 


1.4 


0.0410 


0.0414 


0.0841 


-0.0761 


0.0959 


1.8 


0.0184 


0.0183 


0.0574 


-0.464 


0.0791 


2.2 


0.00823 


0.00822 


0.0370 


-1.231 


0.0939 


2.6 


0.00368 


0.00368 


0.0216 


-2.434 


0.189 


3.0 


0.00165 


0.00165 


0.0116 


-3.969 


0.521 


3.4 


0.000734 


0.000735 


0.00636 


-5.594 


1.402 


3.8 


0„000333 


0.000327 


0,00416 


-7.024 


3.246 


4.2 


0.000153 


0.000145 


0.00350 


-8.028 


6.407 


4.6 


0.0000667 


0.0000643 


0.00334 


-8.482 


11.003 


5.0 


0.0000300 


0.0000283 


0.00318 


-8.387 


16.817 
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In comparing these various solutions, it is noted that the approximate 



transform solution conforms quite favorably with the exact solution using 
up through five terms of the series of equation (4-6) . Truncating beyond 
five terms, the series has a tendency to diverge, the direction and mag- 
nitude of divergence depending on the number of terms beyond five that 
are utilized. Since this series should converge for all values of t, it 
can be concluded that the original derivation of the transform: 



„ r ,. 4 , x( 0 ) x( 0 ) X(0) , 

£ Lx(t) J = + 2 “- + 3 ^ + ... is invalid for large values of t and 

s s s 

is the direct cause of this divergence. For small values of time, t, 



however, all solutions closely approximate the true solution. 



B. EXAMPLE OF A NON-LINEAR EQUATION WITH SQUARE OF DERIVATIVE, 
EQUATION 2 

Considered now is the Froude Equation in its general form: 
x + 2£ (1 + a x) X + x = 0 (4-7) 

This equation demonstrates non-linear damping, but otherwise behaves 
in a rather linear fashion. For a stable system select f = 0 . 1 , a = 0 .5 . 

Then equation (4-7) becomes: 

X + . 2x + . 1(X) 2 + x = 0 (4-8) 

Taking transforms we have: 

s s X(s) - sx(0) - X(0) + .2sX(s) - ,2x(0) + .ls 3 X 2 (s) - .2x(0)s s X(s) 

+ . l[x(0)] 2 s + X(s) = 0. (4-9) 

Collecting terms and letting x(0) = 1 and x (0) = 0, 

. ls 3 X s (s) + [.8s 2 + .2s + 1 ] X(s) - (.9s + .2) = 0 (4-10) 
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Now solve for X(s) using the quadratic formula and arrange in a 



suitable form: 



x(s) . L3^l±J1 } (4-i D 

Expanding the term within the radical in a binomial series and after 
collecting terms and rearranging, X(s) can be expressed in a series. 



w q \ = ( ._9 s + . 2 ) 

1 ’ (.8s 2 + .2s + 1) 



,ls a (.9s + ,Z) S 



,02s s (.9s + ,2) ; 



(.8s^+.2s + l) J (. 8 s 2 + ,2s + 1) e 



(4-12 



The inverse Laplace Transform, x(t), is now obtained using the 
method depicted in Appendix A. The solution, compared with the equa- 
tion solution using a Runge-Kutta integration method, is displayed in 
Table IV'- 0 . 

The Runge-Kutta (true) solution and the approximate solution with 
values corresponding to those displayed in Table IV-2 are depicted in 
graphical form in Figure 1. The true solution is a damped sinusoid. The 
approximate solution in this case is fairly accurate for small values of 
t (up to 3 seconds). This conforms with the analysis previously given in 
Section IV-A. 

Introducing appropriate initial conditions at periodic intervals of 
approximately three seconds , several piecewise but continuing solutions 
were obtained. This method provided a continuous solution but one that 
could be obtained over short periods of time. Figure 1 shows that a very 
accurate solution could be obtained in this manner. 
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TABLE IV-2 



SOLUTION TO EQUATION 2 



TIME 


RUNGE-KUTTA 

SOLUTION 


TRANS FORM 
SOLUTION 


0.0 


1.000 


1.000 


0.4 


0.922 


0.899 


00 

• 

o 


0.709 


0.735 


1.2 


0.397 


0.466 


1.6 


0.038 


0.121 


2.0 


-0.311 


-0.253 


2.4 


-0.596 


-0.595 


2.8 


-0.774 


-0.769 


3.2 


-0.823 


-0.807 


3.6 


-0.743 


-0.707 


4.0 


-0.558 


-0.493 


4.4 


-0.308 


-0.207 


4.8 


-0.037 


0.054 


5.2 


0.212 


0.205 


5.6 


0.405 


0.400 


6.0 


0.520 


0.520 


6.4 


0.546 


0.557 


6.8 


0.487 


0,515 


7.2 


0.358 


0.407 


7.6 


0.182 


0.249 


8.0 


-0.012 


0.061 
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PIECEWISE SOLUTION OF EQUATION 2 




m 
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TIME (SEC) FIG. 



CONTINUOUS SOLUTION OF EQUATION 2 




(*)x 
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TIME (SEC) FIG. 



In contrast to this method Figure 2 depicts a true and approximate 
solution, the latter being obtained in a continuous fashion from zero to 
eight seconds. It is noted that prior to three seconds of time the solution 
is fairly accurate but subsequent to this time begins to diverge and 
becomes slightly out of phase. The divergence is noted to increase with 
time . 

C. EXAMPLE OF THE PRODUCT OF THE VARIABLE AND ITS DERIVATIVE, 
EQUATION 3 

We next consider the van der Pol Equation which describes various 
physical situations , the more noted representing an electronic oscillator. 
The general expression for this equation is: 

Jc - a (1 -x 3 )& + kx = U 

Where U represents the input and a and k are contants. Consider the 
system with zero input and with k =1 , then the equation becomes: 

M - a & + a x 3 & + x = 0 (4-12) 

Apply derived transforms. 

s 3 X(s) - sx(0) - x(0) -as X(s) + a x(0) + a s 3 X 3 (s) - a s 3 x(0)X 3 (s) 

+ X(s) = 0 (4-13) 

After collecting terms the expression becomes: 

as 3 X 3 (s) - as 3 x(0)X 3 s + [s 3 - as +l]X(s)-[sx(0)- ax(0) +x(0)]=0 

(4-14) 

Solving this cubic equation for X(s) is no easy task; however, if we 
assume that the initial conditions for the higher derivative terms are 
negligible and equally zero then: 
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